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ASCA OBSERVATIONS OF GROUPS AT RADII OF LOW OVERDENSITY. 
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ABSTRACT 

Through a 3d-modeling of ASCA observations, we performed a spatially resolved X-ray spectroscopic study, 
extending to radii exceeding 150 kpc, for a sample of 9 groups of galaxies. Combined with published ROSAT 
results, we conclude that these systems generally exhibit a strong temperature decline at outer radii. In our best 
case, NGC3268, this corresponds to a flattening of the entropy profile at a level of ~ 400 keV cm 2 . This value is 
high compared both to the observed entropy floor of ~ 100 keV cm 2 and to the expected value from gravitational 
heating. We suggest that the observed entropy profile in most groups at densities exceeding 500 times the critical 
is purely driven by non-gravitational heating processes. After comparison with a larger sample of groups and 
clusters, we conclude that there is a variation in the level of non-gravitational heating between ^100 keV cm 2 
and ~ 400 keV cm 2 within every system. Using models of cluster formation as a reference frame, we established 
that the accreted gas reaches an entropy level of 400 keV cm 2 by redshift 2.0-2.5, while such high entropies 
where not present at redshifts higher than 2.8-3.5, favoring nearly instantaneous preheating. Adopting galactic 
winds as a source of preheating, and scaling the released energy by the observed metal abundance, the variation 
in the preheating could be ascribed mostly to variation in the typical overdensity of the energy injection, ~ 30 for 
an entropy floor (100 keV cm 2 ) and to ~ 5 for an entropy of 400 keV cm 2 . 

Subject headings: Abundances — galaxies: clusters: general — galaxies: evolution — intergalactic medium — 
stars: supernovae — X-rays: galaxies 



1. INTRODUCTION 

Comparative studies of the scaling relations in clusters of 
galaxies reveal strong deviations of the observed relations from 
predictions based on self-similar collapse (Evrard & Henry 
1991; Bower 1996; Loewenstein 2000; Finoguenov, Reiprich, 
Bohringer 2001). These deviations are best characterized by an 
entropy floor in the X-ray emitting gas. While negligible com- 
pared to the entropy due to the accretion shock in large clus- 
ters, such preheating 1 leads to very extended, low-density gas 
distributions on the scales of groups, that cause a steepening of 
most scaling relations (Ponman, Cannon, Navarro 1999, here- 
after PCN). When the temperature of the preheated gas reaches 
the level comparable to the virial temperature, simulations pre- 
dict adiabatic accretion of the gas (e.g. Tozzi, Scharf, Norman 
2000). This also implies strong temperature gradients, propor- 
tional to the gas density to the 2/3 power (7 = 5/3). As noted 
by Loewenstein (2000), the observation of this phenomena is a 
critical test for the importance of preheating, as well as other 
characteristics of cluster collapse. So far no clear examples of 
systems with an adiabatic gas distribution have been found. 

Throughout this Paper we will define the entropy as S = 
kT e /n 2 J 3 h-^ 3 in keV cm-, following PCN. The observational 
data is presented assuming H Q = 50 km s" 1 Mpc" 1 . 

2. OBSERVATIONS 

We select our sample of groups by cross-correlating ASCA 
(Tanaka, Inoue & Holt 1994) and ROSAT (Truemper 1983) 
archival data on elliptical galaxies with the group catalogs by 

'We define preheating as a change in the initial adiabat for the gas before its accretion into the potentials of groups and clusters. The plausible sources of preheating 
include galactic winds and AGN activity, as well as gravitational shocks associated with the formation of large scale structure. 



Huchra & Geller (1982) and Hickson (1982). From an initial 
sample of 23 groups, we omit 12 because of insufficient cov- 
erage by ASCA SIS (Burke et al. 1991) or availability of only 
ROSAT HRI data, that after a critical check, does not allow 
the study of the surface brightness profile at the radii of inter- 
est. Specifically, our criterium was source extent exceeding 150 
kpc. In addition, we exclude HCG92 and NGC3923, because 
the 'group' component of these two systems does not appear 
clearly in the ASCA data (for HCG92, see also ROSAT HRI 
analysis of Pietsch et al. 1997). 

All observations were screened using FTOOLS version 4.2 
with standard screening criteria. An advantage of ASCA data 
is its low background, compared to Chandra and XMM, which 
is essential for studies of low surface brightness emission. An 
obvious drawback is the complexity of the broad ASCA PSF 
Detailed treatment of the PSF effects is therefore a critical part 
of our analysis, where we follow the approach described in 
Finoguenov et al. (1999), also including the geometrical pro- 
jection of the three-dimensional distribution of X-ray emitting 
gas. Compared to clusters of galaxies, however, the temper- 
ature structure of groups is defined by the Fe L-shell peak, 
thus decreasing the importance of energy dependent effects of 
ASCA PSF on temperature determination, as they affect pri- 
marily the continuum. In addition, for energies below 2 keV, 
the energy dependence of the ASCA PSF is negligible. Our 
minimization routines are based on the x 2 criterion. No en- 
ergy binning is done, but an error calculation is introduced as 
in Churazov et al. (1996) to properly account for small number 
statistics. Model fits to the ROSAT surface brightness profiles 
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of these groups are used as the input to the ASCA data model- 
ing. The details of our minimization procedure for ASCA spec- 
tral analysis are described in Finoguenov and Ponman (1999). 
We adapted the XSPEC analysis package to perform the ac- 
tual fitting and error estimation. The spatially resolved spectral 
characteristics are quoted as the best fit solution, plus an esti- 
mate of the 90 % confidence area for possible parameter vari- 
ation, based on the regularization technique (Press et al. 1992; 
Finoguenov and Ponman 1999). To study the systematic errors 
related to the spatially resolved spectroscopic analysis of the 
ASCA data, we follow the approach described in Finoguenov, 
David and Ponman (2000; hereafter FDP). For all ROSAT 
imaging analysis, we use the software described in Snowden 
et al. (1994) and references therein. 

We use the MEKAL plasma code (Mewe et al. 1985, Mewe 
and Kaastra 1995, Liedahl et al. 1995) in all of our spectral anal- 
ysis. Spectral fitting was twofold. First we use the energy range 
0.8-7 keV and determine the intensity of the hard component, 
arising from the unresolved point sources in the central galaxy, 
fixing its spectral shape to a bremsstrahlung temperature of 6.5 
keV (Finoguenov & Jones 2000). This is mostly important for 
heavy element abundance determination (Matsushita 1998) and 
is a specific of the analysis of X-ray spectra of early-type galax- 
ies, where the gas-to-light ratio is an order of magnitude lower 
than in clusters of galaxies. Later we perform the analysis in the 
0.8-3 keV band, where the group's emission is clearly seen over 
the background. The exception is NGC3268, which is much 
hotter (2.5 keV) than other groups. For it we use the 0.8-7 keV 
energy band. Below 1.5 keV, we add a 10% systematic error at 
the 68% confidence level to the spectra to account for system- 
atics in spectral modeling. 

The groups, centered on NGC2300, IC4296, NGC5129, 
NGC6329 and NGC2563 have in the ASCA fields of view 
point-like background sources, which are brighter than the 
'blank-sky' threshold chosen for ASCA background files and 
thus should be removed individually. In these cases we extract 
the spectra of point sources using both the ROSAT/PSPC and 
the ASCA/SIS, subtract the background, analyze the spectra 
and then, at all energies of the ASCA response matrix, esti- 
mate the contribution from these point sources to the regions of 
diffuse emission selected for the analysis. 

Basic characteristics of the sample are listed in Tabjl]. Col- 
umn (1) identifies the system, (2) is the adopted luminosity dis- 
tance, (3) the corresponding scale length. Columns (4-5) give 
the results of the surface brightness fitting outside the central 
region using a j3 model. Column (6) gives the best fit clus- 
ter emission-weighted temperature, (7) the corresponding 90% 
error, and (8) an estimate of the virial radius of the system 
(r m = 1.237^*50 Mpc; Evrard, Metzler, Navarro 1996) us- 
ing kT e . Estimates of the virial radii using X-ray observations 
by Finoguenov, Reiprich and Btihringer (2001) show that the 
actual virial radii are 20% lower than the estimation of Evrard 
et al. (1996). Columns (9-12) list the luminosity in the 0.5-2 
keV band and a fractional contribution from the line emission 
for the system in total and the central galaxy (excluding the hard 
component). The deprojected luminosity is calculated using the 
profiles of temperature, element abundance and normalization, 
obtained in our analysis of ASCA SIS data. The luminosity is 
not corrected for differences in normalization between ASCA 
SIS and other experiments. Since line emission in our energy 
band is dominated by iron, in calculating the fractional contri- 
bution of the line emission, we use the derived iron abundance 
profile and assume a solar abundance ratio for other elements. 



We attribute to the central galaxy the X-ray emission of temper- 
ature matching the stellar velocity dispersion, which is typically 
one-two central bins in our analysis. 

3. RADIAL PROFILES OF TEMPERATURE AND ENTROPY 

Temperature and entropy profiles are presented in Figs. |l] and 
|2| Heavy element abundances are of secondary interest to this 
paper and are discussed in the appendix. The data on temper- 
ature, entropy, Fe and Si abundance is given in Tab.|| The 
general behavior of the temperature profiles is defined by kine- 
matics of the stellar mass loss at the center (Matsushita 2001), 
characterized by the velocity dispersion of stars in the cen- 
tral galaxy. Beyond radii of 150 kpc, the galaxy contribu- 
tion is small and the emission is dominated by the intragroup 
gas, which exhibits a gradual decrease in temperature, extend- 
ing to the observational limits. In our sample, the best exam- 
ples of such temperature behavior are the NGC3268, NGC5 129 
and IC4296 groups. A central drop in temperature is not seen 
in NGC3268, suggesting that the emission is dominated by 
the ICM gas at radii as small as 50 kpc. ROSAT observa- 
tions, extending to larger radii, show that the temperature also 
drops in NGC2563, NGC507 (Buote 2000), and possibly in 
NGC2300 (Davis et al. 1996). A decrease of the temperature 
profile at radii exceeding 150 kpc is commonly seen in groups 
(Trinchieri et al. 1997; Finoguenov et al. 1999; Finoguenov & 
Ponman 1999, Buote 2000, Finoguenov & Jones 2000). Total 
gravitational mass estimates for the five best observed systems 
(NGC3268, NGC6329, NGC4325, IC4296 and NGC5129) are 
given in Finoguenov, Reiprich, Bohringer (2001). 

In deriving the temperature profiles, we were particularly in- 
terested in revealing any isentropic zones in groups. As one can 
see from Fig.|l], where temperature profiles are shown, several 
groups show sharp temperature declines at large radii, typically 
exceeding 200 kpc. Thus the detection of isentropic behavior is 
limited by the spatial extent of the data. 



Isentropic outskirts in the NGC3268 group 

The X-ray emission from the Antlia (NGC3268) group is ex- 
tended, filling the ROSAT PSPC field of view. The determi- 
nation of the surface brightness with the ROSAT PSPC has a 
large uncertainty, due to the source extent and uncertain level 
of the ROSAT background, caused by the short exposure time 
(6 ksec). Inclusion of ASCA data for the surface brightness 
analysis constrains the (3 parameter, because ASCA has a lower 
background at the radii of interest, due to additional pointings 
that cover the group outskirts. With four pointings, about 30 
ksec each, the ASCA data cover the central 35' of the cluster. 
The results of the analysis are listed in Tab.|l]. 

Given the flat gas density profile, projection effects are im- 
portant. ASCA SIS projected temperature is ~ 2 keV, in agree- 
ment with the GIS results of Nakazawa et al. (2000), while the 
deprojected temperature at the center is ~ 2.7 keV (Figjlj). We 
note that the temperature drop from 0.2 to 0.6 Mpc in the Antlia 
cluster is detected in two off-set ASCA pointings, that cover the 
south-west part of the cluster. From the constant entropy level 
seen at outer radii in the Antlia cluster (Fig.||), we conclude that 
the temperature drop is primarily due to the cluster emission 
and not to a contribution from NGC3258, located south-west 
of NGC3268. Our arguments are based on a factor of 3 lower 
entropy, typical of the group centers (Lloyd-Davies et al. 2000 
and this work), compared to the adiabat of the gas at the place of 
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FIG. 1 . — Deprojected temperature profiles, derived using a single temperature model to fit the data. The solid lines correspond to the best-fit with filled circles 
indicating the spatial binning used in the analysis. Dark and light shaded zones around the best fit curves denote the 68 and 90 per cent confidence areas. Black 
crosses in the NGC507 and NGC2563 panels show the results of a similar deprojection analysis applied to ROSAT PSPC data by Buote (2000). A black cross in the 
NGC2300 panel indicates ROSAT PSPC results for NGC2300 outskirts (Davis et al. 1996) and similarly for NGC7619 from Trinchieri, Fabiano, Kim (1997). The 
confidence level for the cited results corresponds to 68%. 
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Deprojected luminosity in units of 10 41 ergs/s, in the 0.5-2 keV band 

also referred to as the Antlia cluster or NGC3258, the second brightest galaxy 

fraction of line emission, based on the iron abundance and assuming the solar abundance ratio 
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FIG. 2. — Entropy profiles in units of keV cm 2 , derived using a single temperature model to fit the data. The error bars are shown at 90% confidence level. 
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FIG. 3. — Left panel. Entropy profiles scaled by the total mass in units of 10 13 Mq, measured at an overdensity of 500 and plotted against a scaling radius measured 
as a fraction of r 5m . Black solid lines denote systems with kT< 1.5 keV (NGC5129, IC4296, NGC4325, HCG62, NGC5044, HCG51, NGC6329), black dotted lines 
1.5 keV <kT< 3 keV (MKW4, NGC3268, A262, MKW4S, A539, AWM4, MKW9, A2634), gray dashed 3 keV <kT< 5 keV (A1060, A4038, 2A0335, HCG94, 
A2052, AWM7, A2657, MKW3S, A2063, Hydra A, A3112, A4059, A496, A2199), gray dotted lines 5 keV <kT< 7.5 keV (A3558, A119, A1651, A3571) and 
gray solid lines 7.5 keV <kT< 11 keV (A401, A3266, A2029). Right panel. Scaled entropy at the overdensity of 500 vs the total gravitational mass of the system, 
measured within an overdensity of 500. The short-dashed line indicates the prediction from shock heating, corrected for the mean formation redshift, normalized to fit 
the intermediate-mass clusters. Thinner dashed lines indicate a scatter induced by the expected variation in the formation redshifts on the 75% level. The long-dashed 
line indicates the effect of the preheating value of 500 keV cm 2 . 



NGC3258. In addition, enhancement in the surface brightness 
would be seen in the ROSAT map, as well as more than just 
one galaxy (NGC3258) ought to be associated with the puta- 
tive subgroup. Finally a coincidental match of the temperature 
close to the center of the Antlia cluster and at larger distances, 
reaching 600 kpc, would require fine tuning. For another group, 
AS 1101, first XMM-Newton results reveal a temperature pro- 
file in remarkable correspondence to NGC3268 (Kaastra et al. 
2001). 

From radii of 150 kpc to 600 kpc, the entropy profile of the 
Antlia cluster is flat at the level of 430 keV cm 2 . In many groups 
the temperature drop is very sharp, corresponding to the level- 
ing off of the entropy at a value of 350 ± 50 keV cm 2 (see Fig.[|). 
While the plateau in the entropy profile is predicted by the pre- 
heating scenario, the observed entropy level is 3 times higher 
than the measurement of PCN. In the following we consider 
the reasons for this phenomenon. 

4. A CEILING IN THE PREHEATING VALUE FOR ENTROPY 

We compare the observed entropy level at the group's out- 
skirts to that predicted from scaling relations, choosing an over- 
density (75/ pent - 1) of 500, A500, where ~p is the mean density 
of the cluster integrated out to r& . In Finoguenov, Reiprich and 
Bohringer (2001) we estimated radii for A500 (hereafter rsoo), 
the polytropic index and the total mass M 5 oo, while in Sander- 
son et al. (2002) we fit a /3-model to the observed ASCA de- 
projected density. The employed /3-model fits to surface bright- 
ness profiles and polytropic fits to the temperature profiles are 
only valid outside the cores of the systems, typically larger than 
O.lrsoo- The resulting entropy profiles for samples of groups 
and clusters are plotted in FigB vs the fraction of r 5 oo, scaled by 



A^oOj . the total enclosed mass measured within A500 in units of 
10 13 Mq. Such scaling is equivalent to scaling by virial temper- 
ature and is introduced to remove the dependence of the shock 
strength on the cluster mass. We use the mass instead of the 
measured temperature to exclude the effect of adiabatic com- 
pression. Entropy profiles produced by the accretion shock ex- 
hibit the strongest rise with radius, as seen for rich clusters (also 
David et al. 1996). Entropy profiles of groups are much shal- 
lower (solid and dotted lines in Fig.||) and, especially at smaller 
radii, are higher than found for clusters. 

The flatter entropy profiles found for groups are in good 
agreement with simulations that introduce an entropy floor 
(Borgani et al. 2001). However, at low overdensities groups 
still deviate from clusters. These deviations are opposite to the 
prediction of structure formation, since groups form earlier than 
clusters, resulting in lower entropy. These deviations are larger 
than expected due to the entropy floor (PCN) and the effect of 
gas cooling (Voit & Bryan 2001). 

4. 1 . Effects of gravitational heating and substructure 

Clusters do not appear as uniform entities. As Fig.|| (left) 
shows, at 0.5r5oo the entropy profiles of clusters start to diverge, 
so the scatter in the entropy among clusters at a given over- 
density becomes similar to the differences between clusters and 
groups. A critical radius at which such deviations occur corre- 
sponds to rsoo, matching an earlier suggestion of Evrard et al. 
(1996) that X-ray determinations of cluster mass should gen- 
erally be trusted to such radii only. We produce a cut in the 
entropy profiles at 7-500 and in Fig.g (right) plot the resulting en- 
tropy values vs the gravitational mass of the system, measured 
at overdensity of 500. Preheating at the level discussed in this 
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paper (shown as a long-dashed line) can only affect systems 
with mass (M500) below 10 14 M©, and is successful in describ- 
ing the entropy behavior in groups. 

In Fig.|] (right) we also present the rescaled theoretical pre- 
diction for shock heating (Eke, Navarro, Frenk 1998; PCN), 
with details presented in the appendix. As is seen from the Fig- 
ure, clusters exhibit a large spread around the theoretical predic- 
tion, compared to the effect of variation in the cluster formation 
epochs. 

In this paragraph we consider the possible reasons for the 
large spread in the entropy among clusters. Incomplete merger 
shock heating predicts a sharp temperature and density changes, 
arising from a final time of shock propagation through the clus- 
ter. Since Chandra cluster images often show surface brightness 
edges (Markevitch 2002), the temperature behavior generally 
suggests cold fronts rather than shocks (Markevitch et al. 2000; 
Vikhlinin, Markevitch & Murray 2001). Another possibility 
is that merger shocks start in the cluster core and decrease in 
strength with radius (due to the flat distribution of the gas den- 
sity). Finally, survival of the gas associated with the infalling 
subcluster results in a decreased entropy (e.g. Borgani et al. 
2002). Since formation of large clusters is still in progress, 
compared to poor clusters (Kauffmann 1995), hotter clusters 
tend to have lower entropy, resulting in a larger range around 
the entropy prediction which includes the effect of cluster for- 
mation epoch. In addition, an observed scatter in the entropy 
could result from our simplified approximations to both tem- 
perature and density profiles. At least in one example of the 
hot cluster (although not from our sample), A2163, where a de- 
creasing temperature profile from ASCA is confirmed by Chan- 
dra and XMM-Newton (Markevitch, Vikhlinin 2001; Pratt, Ar- 
naud, Aghanim 200 1 ), a steepening of the density profile is also 
suggested at large radii (Vikhlinin, Forman, Jones 1999a). If 
such behavior is present in other systems, the resulting entropy 
would be underestimated. 

4.2. Leaning over the gravitational threshold 

At radii beyond ~ O.lrsoo, entropy profiles in clusters may 
be the result of purely gravitational heating. Preheating will in- 
crease the scatter in the entropy profiles between clusters and 
groups. One can see from Fig.|| (left) that the scatter between 
cluster and group profiles, compared to the scatter within the 
cluster profiles, decreases with increasing radius. For example 
at half rsoo, the range among all the systems is a factor of 10, 
while it is a factor of 2 for clusters. At r$oo the scatter is a 
factor of 3 for all systems and a factor of 2 for clusters, while 
at yet larger radii the scatter does not change when groups are 
excluded. This indicates that r^oo is a limiting radius where 
the scatter arising from the non-gravitational effects can still be 
seen above the scatter resulting from the gravitational heating. 
This corresponds to a typical mass for the groups in our sample 
of around 3 x 10 13 M Q . 

Let's compare the entropy level measured for radii of equal 
enclosed mass for groups and clusters. Observations show that 
there is no significant trend in gas fraction with system mass 
at radii defined this way (Sanderson et al. 2003), which is im- 
portant in estimating the entropy produced by a given shock 
strength, as explained in § |l.6| . Consideration of the entropy at 
similar enclosed mass is fruitful because the implied strength 
of the accretion shock is thought to be similar, unless the epoch 
of cluster formation plays a dominant role. We choose a mass 

2 In this approach we are limited to samples analyzed by Finoguenov (FDP, 
A119, A3571, A401, A3266. 



of 3 x 10 13 M Q for this comparison, which is still within the 
analyzed radii for the groups, while outside the 150 kpc radii 
for clusters in our sample (and thus is free from complications 
due to BCG potential, cooling flow and BCG AGN heating). 

A comparison of the entropy for this enclosed mass is pre- 
sented in Fig.|](left panel) for each 2 of the clusters and groups 
(38 systems in total) shown in FigjJ. An entropy level of 350 
keV cm 2 is typical for systems with temperatures below 3 keV. 
A dramatic change in the behavior of entropy vs temperature 
occurs above 3 keV, exhibiting a trend towards lower entropy 
for hotter systems. The dashed line in the left panel of FigJ^ 
shows the entropy attributed to shock heating, calculated for 
the mean formation redshift (z/) of groups and clusters (the 
solid line on the right panel of FigQ) and up-scaled by a fac- 
tor of two (it is off the graph when the normalization of PCN is 
adopted). It is clearly seen that the observed trend is different 
from the entropy behavior predicted from shock heating. 

If we assume a similar history of gas preheating in groups 
and clusters, then the trend in Fig.|] (left) requires a growth in 
the gas entropy with (Hubble) time, with the entropy increasing 
from cluster cores (zj ~ 2.7) to the outskirts of groups (zj ~ 0.). 
The entropy in poor clusters (kT^ 3 keV; z/ ~ 2.0 ) reaches 
a level, comparable to that found in groups. This picture im- 
plies both prompt preheating at high-z and a lack of preheating 
at low-z. Thus we suggest that the preheating level, observed 
through the entropy, will not increase in other observational 
tests and corresponding therefore to the entropy ceiling. 

An entropy floor has been suggested based on the observation 
of excess entropy in cores of groups compared to the gravita- 
tional heating (PCN). Introduction of the term "ceiling" would 
represent the upper limit on the preheating, which is a charac- 
teristic of any realistic model of the preheating, based on the 
limited energetics provided by putative sources of preheating, 
such as SNe. Term entropy ceiling is also appropriate for sce- 
narios ignoring the preheating, where one can reformulate a key 
assumption of Voit & Bryan (2001), to read 'a ceiling in the pre- 
heating value for entropy is much lower than the entropy floor' . 
Since we prove otherwise, the scenario of Voit & Bryan (2001) 
in its original formulation is dismissed. Below, we present a 
more detailed consideration of the establishment of the entropy 
ceiling, after first noting other supporting observations. 

The entropy level for the gas accreting onto the outskirts of 
clusters can be estimated from the entropy level of filaments 
connecting the clusters. Durret et al. (1998) found that a galaxy 
filament close to A85 exhibits an elongated X-ray emission in 
the 0.4-2 keV band. From the reported flux of 2.7 x 10~ 13 ergs 
cm" 2 s" 1 , estimating the flux collecting area as 60' x 10', and 
taking the short axis for the third dimension, we obtain a den- 
sity of the X-ray emitting electrons of ~ 3 x lO" 5 /!^ 2 cm" 3 . For 
emission to be detected in the 0.4-2 keV band, the gas tem- 
perature should exceed 0.3 keV, corresponding to an entropy 
of > 310 keV cm 2 . Briel and Henry (1995) used the ROSAT 
all sky survey data to combine the measurements of regions 
connecting physically close clusters of galaxies, identified in 
Huchra et al. (1990). They obtain a 2a limit on the emission 
of 3.9 x 10~ 16 ergs cm" 2 s . With the assumption of a 0.5 keV 
temperature for the gas in filaments, justified by observation 
of metals in clusters of galaxies (Finoguenov, Arnaud, David 
2001; hereafter FAD), Briel and Henry (1995) derived a limit 



and this paper). This excludes a few hot clusters, introduced in FigM A3558, 
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FIG. 4. — Left panel. Entropy of the system derived at equal enclosed mass (3 X 10 13 Mq) plotted against the system temperature. Dashed line represents a predic- 
tion from shock heating artificially scaled up (by a factor of 2) to reproduce the group measurements. Two dotted lines separate the sample into three groups. Systems 
with entropy exceeding 300 keV cm 2 are assumed to achieve the highest preheating level, the systems with entropy below 200 keV cm 2 are counted as forming before 
the high entropy level is achieved and the rest of the systems are considered unconstrained. Right panel. Virialization redshift of the first 3 X 10 13 MQin clusters vs the 
observed temperature. Solid-line curve shows the location of 50% of the systems from numerical simulations of structure formation (Lacey & Cole 1993), calculated 
using the formulae in Balogh et al. (1999) for f2 m , = 0.3. Short dashed-line curves show the location of formation redshift for 10% and 90% of systems. Shaded 
zones on both panels indicate whether or not the high level of preheating has been achieved. Light shadowing indicates that 75% of groups has been formed after and 
two dark shadows indicate that 25% of cluster cores form before and another 25% form after the high level of preheating has been achieved. Two long-dashed lines 
show a variation in the redshift consistent with the data, while two dotted lines indicate the limits where the model prediction becomes strongly inconsistent with the 
data. 



on the gas density 3 of 7.4 x 10 _5 /i 1//2 cm" 3 , which for h = 0.5 
translates into a lower limit on the entropy of 360 keV cm 2 . 
Thus we conclude that also in rich clusters, the entropy of the 
accreting gas is higher than the entropy floor. 

Observations of a newly discovered class of objects, fossil 
groups or OLEGs (Ponman et al. 1994; Vikhlinin et al. 1999b) 
also show a high level of preheating (entropy is 350 keV cm 2 at 
radii of the enclosed mass of 3 x 10 13 M , Vikhlinin, A. 2001, 
private communication). 

Interpretation of the M gas - T relation for clusters of galax- 
ies also requires a higher preheating level of ~ 400 keV cm 2 
(Babul et al. 2002, McCarthy et al. 2002). The issue of the 
level of preheating in the theoretical models is still unsettled, 
as e.g. model of Dos Santos and Dore (2002) does not require 
such a high entropy. 

4.3. Dating the preheating using a prescription for cluster 
formation 

With recent advances in high-resolution numerical simula- 
tions, cosmic microwave background observations, the SN la 
Cosmology project, and surveys of clusters of galaxies, stud- 
ies of the evolution of the Universe have achieved a convincing 
picture with a large degree of detail, involving dark energy and 
dark matter. These constrain the geometry of the Universe, its 
dynamics and the large-scale structure of the dark matter distri- 
bution. Thus, from the difference in the formation epochs for 
the outskirts of groups compared to cluster cores, we can infer 
the redshift over which the entropy ceiling has been maintained. 

3 assuming an iron abundance of 0.3 solar 



In order to estimate the redshift when most of the preheat- 
ing occurs, we show in FigJ^ (right) the analytical and numer- 
ical results for the formation of clusters (Lacey & Cole 1993), 
following the formulae presented in Balogh et al. (1999) for 
^m,o = 0.3. With a solid line we indicate the mean of the cluster 
formation redshift distribution and with dashed lines the corre- 
sponding 10% and 90% values. To represent the measurements 
on the left panel, we divide the systems into three groups ac- 
cording to the entropy as low, intermediate and high, < 200, 
200-300, > 300 keV cm 2 , respectively. We attribute the gas 
with entropy lower than 200 keV cm 2 to that accreted before the 
establishment of the entropy ceiling, while the gas with entropy 
higher than 300 keV cm 2 is considered to achieve the entropy 
ceiling before accretion. Analysis of the left panel reveals that 
at least 75% of the low temperature systems have achieved a 
high level of entropy and therefore have formed later than the 
entropy ceiling was established. This is shown by the light-grey 
shaded zone in the right panel. The dark-shaded zone at higher 
temperatures indicates the presence of both high (at least 25% 
of clusters) and low (at least 25% of clusters) entropy systems. 
Low entropy systems are treated as formed before the entropy 
ceiling was established. 

As is shown in Fig.|] (right), at redshifts greater than 2.8- 
3.5 the accreted gas was characterized by low entropy, while 
starting at redshifts 2-2.5, the accreted gas is on a high adiabat. 
Duration of this transition period is very short, 0.3- 1.4 Gyr. 
The limits on the epoch for feedback processes, based on the 
entropy measurements, crucially depend on the way such feed- 
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FIG. 5. — Search for a correlation of the entropy measured at radii of equal enclosed mass (3 X 10 13 Mq) vs Fe (left) and Si (right) abundance. Groups are shown 
with filled black circles, cool clusters (2 keV <kT < 3.5 keV) with filled grey circles and hot clusters with open circles. 



back occurs. However, a lower limit of Zf > 2 on the feedback 
redshift should hold in any scenario for external (pre) heating. 
Yet, the requirement of high redshift for the feedback processes 
responsible for the preheating further supports the dominant 
role of bulges and is an important input to modeling galaxy 
formation in high-density environments. This conclusion is in 
agreement with observational evidence of no star-formation in 
early-type galaxies since a redshift of 2 (Ellis et al. 1997; Ben- 
der et al 1998). 

4.4. Relating enrichment to the preheating 

In this section we attribute the establishment of the preheat- 
ing level for the gas to action of the galactic winds. Within this 
scenario, the metals found in the hot intracluster medium can 
serve as an additional constraint on the preheating, as they are 
directly related to the amount of energy released into the ICM 
(FAD). 

In the previous section we concluded that a high preheat- 
ing level corresponding to an entropy of about 350 keV cm 2 
is reached at large radii in all systems. On the other hand, 
the observed level of heavy element enrichment varies between 
groups and clusters leading to the suggestion that galaxy winds 
deposit more energy in clusters (FAD). Moreover, such enrich- 
ment was found to be constant over the cluster (FDP), while a 
simple interpretation of the results of the previous section leads 
to the expectation of an abundance profile rising with radius. To 
better illustrate this point, we show in Fig.^J that the entropy at 
a radius of equal enclosed mass does not correlate with either 
the Fe or Si abundance. 

Apparently, at least part of the enrichment is not associated 
with the preheating {e.g. SN la enrichment occurs at later 
times, which determines the Fe abundance in clusters at radii 
corresponding to Fig.||). Furthermore, Chandra and XMM- 
Newton observations do not significantly change the picture of 
the distribution of SN II products (David et al. 2001 ; Matsushita 



et al. 2002; Molendi & Gastaldello 2001; Tamura et al. 2001). 

We have concluded that preheating of the gas reaching the 
virial radius should be complete by a redshift of 2-3.5. Since 
it takes some time before the preheated gas will be accreted, the 
redshift of the feedback is somewhat higher. Let's choose a red- 
shift of 3.5 for the energy injection and calculate the overden- 
sity at which energy injection takes place in order to reproduce 
the observed entropy for radii corresponding to 3 x 10 13 M Q 
enclosed total mass (it will be possible for us to check the va- 
lidity of a particular assumption for the redshift). We account 
only for the effects of star-formation, estimating the amount 
of energy from direct observations of the chemical enrichment, 
as calculated by FAD. (For simplicity we will take an average 
for systems with similar total mass, which is 0.4 keV/particle 
for groups, reaching 1 .2 keV/particle for rich clusters, and an- 
alytically expressed as 0.44 x (kT /keV) 5 keV/particle). In 
Fig |(] (left) we determine the overdensity where this amount 
of heat would lead to the observed entropy level. We assume 
the electron density at critical density of 4.1 x 10~ 7 cm -3 (tak- 
ing n e /n ion = 1.09, m ion = \ .21m p , f haryon = 0.17). For the cho- 
sen redshift of 3.5, the scaling suitable for f2 m . = 0.3 results in 
«7 = 2.1 x 10~ 5 cm~ 3 corresponding to the critical density. There 
is a slight underestimate of the overdensity since we take the 
density at which we eject the heat as the overdensity. We can 
take into account at least galaxies embedded in the preheated 
gas, if we introduce an increase in the overdensity by 20% due 
to a gas-to-bulge ratio of 5. Variation in the assumption for the 
redshift changes both left and right panels of Fig.|| such that 
matching between the temperature of the system at turn-around 
on both panels is maintained for a broad range of redshifts. A 
requirement for the minimum derived overdensity to be higher 
than linear approximation to the growth of primordial density 
fluctuations (1.69/(1 +Zmjectwn)) rules out redshifts > 4. 

Fig.|| (left) shows that many systems, including the outskirts 
of groups, are characterized by SN enrichment at turn-around, 
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FIG. 6. — Left panel. Overdensity at which galaxy winds have been ejected, calculated from the requirement to match both amount of energy associated with 
metal enrichment as well as the observed entropy level (at radii of 3 X 10 13 Mq enclosed total mass), assuming redshift of star formation of 3.5, and using the scaling 
appropriate for the f2 m o = 0.3. The dramatic difference in the entropy between groups (kT < 2 keV) and clusters, translates into the statement that preheating of 
groups is complete at the turn around radius, while preheating of cluster cores occurred at higher overdensities (yet outside the virial radius). Three dashed lines mark 
overdensities indicating three stages of cluster formation, linear growth of initial density fluctuation (1.69/(1 +z) = 0.38), turn around (4.6) and virialization (150). 
Right panel. Redshift of turn around, derived by taking half of the formation time from FigM Grey line shows a consistency with the assumption in the left panel on 
the redshift of energy injection. Dotted lines indicate the temperature range for systems at which the transition in the derived overdensities occurs. It matches well the 
analytical prediction, shown by solid line. Dashed lines show the expected spread in turn around redshifts, drawn at 75% confidence level. 



while cluster cores are characterized by enrichment occurring 
primarily at higher overdensities. In Fig.|| (right) we present 
a calculation for the turn-around redshift that corresponds to 
the virialization epoch, shown in Fig.Q (right). The turn-around 
time is calculated as half the collapse time, defined in Fig.Q 
(right). This plot allows us to check our assumption for the 
star-formation redshift. For the chosen redshift of 3.5, the gas 
presently accreting onto the outskirts of groups should not nec- 
essary be inside the turn around radius. Thus, there is no con- 
tradiction in our estimate that injection overdensities are less 
than 5. In our approach, the preheating level in groups is de- 
fined by the overdensity of the material destined to collapse 
(0.4 < S < 4.6). Constant entropy behavior at radii of equal 
enclosed mass is therefore a consequence of a narrow range of 
possible overdensities and a small scatter in the metal produc- 
tion. For the cores of clusters with temperatures above 4.5 keV, 
the high overdensity shows that the metal enrichment occurred 
during the advanced stages of cluster core formation. A more 
refined approach, would allow for the expectation that most star 
formation takes place at larger redshifts in the regions which 
now form cluster cores, compared to the outskirts of groups. 
This will reduce the range in derived overdensities at which pre- 
heating is inferred to happen, and in particular it will increase 
the very low overdensities inferred for group outskirts. Another 
effect arises from the importance of SN la in the production of 
metals in groups, generally understood as enrichment at high 
density regions, thus reducing the level of preheating. Both ef- 
fects just mentioned are of the same factor of few, compensating 
each other. Two other uncertainties related to our calculation 
are the energy-to-temperature transformation, as discussed in 
Lloyd-Davies et al. (2001), which we assume to be 3/2kT and 



the energy released per SNe, which we take as 1.5 x 10 51 ergs 
in our calculation. 

In summary, a gradual increase in the entropy of the accreting 
gas in any system could be characterized by a single preheating 
epoch, with the cores of systems being already at an advanced 
stage of collapse. Flat a-element abundance profiles observed 
in clusters of galaxies (FDP) are naturally reproduced, since at 
a derived range of overdensities there is no segregation between 
gas and galaxies. A smoking-gun to define the epoch and level 
of preheating would be an observation of the thermodynamic 
state of baryons in the infall regions of clusters, either in emis- 
sion or absorption. 

4.5. Comparison to existing approaches 

Borgani et al. (2002) have carried out simulations where var- 
ious preheating entropy levels are reproduced. One difference 
between their approach and ours lies in the attitude regarding 
the high density regions. Borgani et al. (2002) rises the entropy 
of all particles in the system to the same level, which requires 
a much higher level of energy injection, as well as preferen- 
tial deposition of metals into high-density regions. The highest 
overdensity where we introduce preheating is 30. Since in our 
approach, we do not add energy to high density regions, we al- 
low star-formation to continue in Ly break galaxies. Still, we 
are able to preheat the bulk of the gas mass, which resides at 
lower overdensities. As a consequence, our approach is dif- 
ferent to that of Borgani et al. (2002) by requiring the dense 
regions to cool, either by forming galaxies outside the virial ra- 
dius or via early cooling flows, as proposed by Voit and Bryan 
(2001). An advantage of our approach lies in lowering the re- 
quirements for the total energy required to achieve the observed 
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preheating level. 

One firm conclusion that can be derived from our work is 
that for star-formation to be responsible for the preheating in 
groups and clusters, the energy ejection should occur at an 
early epoch and preferentially in a low-density environment. 
As star-formation itself usually peaks in high-density regions, 
it might seem that such a requirement could rule out this sce- 
nario. However, examples of local star-formation show that 
galaxies are very good in building chimneys expel most of the 
energy and metals (self-regulating star-formation; Larson & 
Dinerstein 1975). Also if energy is equally deposited at vari- 
ous densities, the low-entropy gas experiences large radiative 
losses and falls back onto the galaxy and is consumed in subse- 
quent star-formation. Since high overdensities surrounding the 
galaxy may not become part of the ICM (Larson & Dinerstein 
1975), it is quite logical to assume that the energy associated 
with metals in the ICM comes from the high-entropy tail of the 
energy ejection (FAD). 

Injecting the energy into overdensities lower than 5 requires 
that the infall of the gas into the group horizon occurs after the 
gas is preheated. Also, when we calculate the amount of metals 
in groups, we only account for the filaments that fall in and not 
those that do not. This explains why it is reasonable to assume 
no losses due to the removal from galaxies by galactic winds of 
the part of metals found in clusters. 

4.6. Alternative sources of preheating 

To estimate the uncertainties in using entropy to constrain 
the feed-back epoch, one can consider the influence of cooling 
flows and AGN heating, as well as segregation between mass 
and gas inside the cluster virial radius. 

The high level of entropy, associated with the entropy ceiling 
is beyond the reach of models that invoke gas cooling (e.g. Voit 
& Bryan 2001) to explain the 100 keV cm 2 entropy floor. We 
have shown that the high entropy level in the groups' outskirts 
is related to recently accreted gas, while lower entropy levels 
would be more typical of an earlier epoch. One effect of gas 
cooling is to reduce of the entropy of the gas, due to entropy 
carried out by emitting photons. However, most of the cooling 
gas is expected to drop out at the same location as the cooling, 
so part of the information originally recorded in the entropy at 
centers of groups and clusters is lost due to the conversion of 
the X-ray emitting gas into stars. Regarding the systematics of 
using only the cores of clusters as probes, note that our systems 
are sampled at a radial range outside the regions, where cooling 
is appreciable and in addition we use only the hotter component 
in our analysis of cluster cores. 

AGNs are proposed as a source of the preheating, as they 
may provide much higher energies per particle, compared to 
SNe (Wu, Fabian & Nulsen 2000). Using the AGN luminosity 
function with a proper scaling by cluster volume and overden- 
sity, and an increase by a factor of 80 in the AGN activity with 
the redshift (Miyaji, Hasinger and Schmidt 2000), the energy 
radiated in X-rays over a Hubble time for all AGN corresponds 
to 0.2 keV/particle. AGN still could play a dominant role if, 
according to Fabian & Iwazawa (1999), most of them are ab- 
sorbed and most of the emission is re-radiated in the infrared. 
Assuming an efficient accretion mode forming 100% of the lo- 
cal volume density of black holes, a corresponding estimate of 
the radiated energy is 10-15 keV/particle. Observationally there 
is a lack of such luminous absorbed AGNs. In particular the 
optical follow-up of the XMM Lockman-Hole field and Chan- 
dra Deep Field South reveal that absorbed AGNs are located at 



low redshifts (Hasinger 2002, Rosati et al. 2002). Although 
the identification of AGNs is not yet complete, the unidentified 
sources are too faint to account for efficient mode of accretion 
growth of the black holes. Ciotti & Ostriker (1999) also argue 
that feedback processes prevent AGN from effective accretion. 

Another problem, related to the suggestion of Fabian (1999) 
on AGN importance is related to the conversion of radiated en- 
ergy to kinetic energy. The most effective mechanism, based 
on radiative pressure (Silk & Rees 1998), requires galactic col- 
umn densities of 10 24 cm 2 , while the bulk of baryons at high 
redshift have column densities less than 10 21 cm 2 (Ellison et al. 
2001). Yamada & Fujita (2001) considered heating by AGN 
jets and concluded that jets are capable of providing about 0.3 
keV/particle at redshifts below 3. They assume that 1% of nor- 
mal galaxies are radio galaxies and limit the redshift to be less 
then 3 using the prediction for the deformation of the cosmic 
microwave background (Sunyaev-Zel'dovich effect). One can 
elaborate this scenario by considering the association of radio- 
loud quasars with elliptical galaxies (e.g. Laor 2000), which 
formally increases the preheating in clusters by up to an or- 
der of magnitude (3 keV/particle). However, Ho (2002) pro- 
posed that the association of radio-loud sources with ellipticals 
is driven by a tendency of local galaxies to have nuclei that 
are both underluminous and radio-loud, in qualitative agree- 
ment with predictions of accretion disk theory. Thus, such an 
association would not hold when ellipticals were star-forming 
and the galactic wind was driven by star-formation. Indeed, as 
acknowledged by Laor (2000), Stern et al. (2000) find no evi- 
dence that the radio-loud fraction depends on optical luminosity 
for -25 > M B > -28 at z = 2-4. 

Another aspect of AGN activity is related to the robustness 
of our approach, where we use the cluster cores. We can con- 
sider Hydra-A as a counter-example, where the AGN activity 
is among the strongest (McNamara et al. 2000) at the present 
time, yet its entropy is amongst lowest (134 ±33 keV cm 2 ). 
As discussed in McNamara et al. (2000), no significant shock 
heating is seen in Hydra-A. The reason is that the typical ki- 
netic energy associated with the observed AGN activity in clus- 
ters centers is on the level of 10 44 - 10 45 ergs/s (Boringer et 
al. 2002), while to substantially change the entropy of the cen- 
tral gas mass of 10 12 M , typical of cluster cores, it requires 
kinetic input from AGN of 10 48 ergs/s. Since the abundance 
of hot clusters is 10~ 8 Mpc" 3 (Reiprich & Bohringer 2002), to 
systematically change the entropy of the cluster cores it would 
require an AGN activity three orders of magnitude higher than 
the 10 -11 Mpc -3 maximal volume density of AGNs with 10 47 
ergs/s, as observed at high redshift. Note that AGNs of such lu- 
minosity are capable of completely ionizing the torus and could 
not be hidden from the observer (Nandra et al. 1995). 

The last possible source of systematics we consider is the 
segregation of mass and gas. In fact the gas fractions are sim- 
ilar between groups and clusters at radii of the same enclosed 
mass and are 0.05h^ 5 (Sanderson et al. 2003) at radii of total 
mass of 3 x 10 13 M Q . Thus the radii of equal enclosed gas mass 
and equal total mass agree in this case. 

Observation of a systematic increase in the gas mass fraction 
with the mass of the system (David et al. 1995; Jones, For- 
man 1999; Sanderson et al. 2002), could be self-consistently 
explained within the preheating scheme (Babul et al. 2002). 
However, if reduction in the gas fraction occurs before the 
shock heating, the entropy resulting from the gravitational heat- 
ing will be higher. This possibility does not change the consid- 
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eration of the entropy at radii of equal enclosed mass, since the 
gas fractions are similar. The gravitational heating fails to re- 
produce the trend in the entropy at radii of equal enclosed mass, 
independent of the normalization, thus disfavoring the leading 
role of the reduction in the gas fraction before the shock in the 
entropy profiles. 

5. CONCLUSIONS 

We present the temperature, entropy and heavy element 
abundance measurements for nine groups of galaxies and com- 
pare their entropy profiles to a larger sample of groups and clus- 
ters of galaxies. We conclude 

• Analysis of groups at radii exceeding 150 kpc generally 
reveals a decrease with radius in the gas temperature. In 
the best-resolved case, NGC3268, this corresponds to a 
flat entropy profile at the level of 400 keV cm 2 . 

• Comparison of the scaled entropy at an overdensity of 
500 between groups and clusters, still reveals deviations, 
that require a high level of preheating (400-500 keV 
cm 2 ). 

• Comparing the entropy in groups and clusters at radii 
with a similar enclosed mass (chosen as 3 x 10 13 M ), 
we conclude that the entropy of the preheated gas does 
not exceed 400 keV cm 2 . Thus, we establish a ceiling for 
the preheating entropy. 

• Comparing our results to the prediction of structure for- 
mation, we find that accretion of the gas preheated to 400 
keV cm 2 did not occur before a redshift of 3.5, but is typ- 
ical of accretion for redshifts below 2. Thus, not only the 
metal enrichment of clusters is associated with bulges, 
but also the timing of the preheating suggests the domi- 
nant role of bulges in the feed-back process. 

• If an increase in the level of preheating is due to secu- 
lar star-formation activity, duration for energy ejection is 
constrained to a narrow time interval (redshifts 2-3.5). 



Distribution of the level of preheating cannot result from 
the progressive release of SN energy from galaxies, but 
can be produced by the nearly instantaneous release of 
SN energy into a density gradient. 

• Assuming an injection (or star-formation) redshift of 3.5 
and scaling the feed-back by the observed metal content 
in groups and clusters of galaxies, (0.4 keV/particle in 
groups, 1.2 keV/particle in clusters), the typical overden- 
sity of gas receiving the preheating must be below 5 for 
outskirts of groups and 5-30 for cluster cores. 

• Using the luminosity function of AGNs, we show that 
AGNs cannot be the dominant source of the preheating. 
A concentration of obscured AGNs at low redshift does 
not support the suggestion of an unseen growth in black- 
hole masses. 

• The "energy crisis", related to the preheating with a lim- 
ited source of energy, such as SNe, is solved in our ap- 
proach by allowing the high-density regions (overdensi- 
ties higher than 50 at redshift of 3.5) to cool out and form 
stars, while we are able to preheat the dominant gas frac- 
tion which resides at low densities. 
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APPENDIX 



A. RESCALING THE ENTROPY FROM GRAVITATIONAL HEATING 

The scaling of the entropy, chosen in Fig. [3] (right) accounts for the increasing strength of the shock with the cluster temperature. 
To scale the theoretical prediction for the shock heating (Eke, Navarro, Frenk 1998; PCN), stated at O.lr,,, we rescale the entropy 
to the radius of the overdensity of 500, using the S ~ r 11 (such entropy scaling on radius is also suggested in Tozzi & Norman 
2000), p m ~ r~ 2A , found in the same simulations (Eke, Navarro, Frenk 1998) and use the M— T relation of Evrard et al. (1996) to 
consistently remove the entropy scaling on temperature. We scale the prediction to h = 0.5, and account for the fact that in the low 
density Universe, the shock occurs at lower overdensity and is even lower when lambda is taken into account (100 compared to 178, 
Pierpaoli, Scott, White 2001), so the shock temperature is lower, but the density is even lower. We obtain 



where M500 is in units of 10 13 M©. 
Compared to the original 



S(/ 5 oo) _.,„„/' fgas 



M, 



500 



189 



0.17 



1+6 

Too 



(i) 



s ^^ 5 {^v){W h ' ikeVcm2 



Contribution of various scaling are following 



S(r 5 oo) = 8.3S(0.1r v ) 



(2) 
(3) 



ft 



= 1.26 



fgas 



'0.5 



^0.06/ V - 17 , 
where the PCN scaling of /T 4 / 3 resulting from p crlt ~ h 1 is balanced by the scaling coming from f gas ~ /T 3 / 2 . 



(4) 



1 + 6 



= 1.47 



1 + 6 



(5) 



178 J V 100 

In Fig.|| (right) we show the prediction for shock heating and include the influence in the redshifts of cluster formation (S ~ 
(1 +Zf)~ l in CDM, while we use a prescription for density and virial overdensity evolution in ACDM). 

B. HEAVY ELEMENT ABUNDANCE DETERMINATION 

All abundances are given relative to the solar photospheric values in Anders & Grevesse (1989). The abundances of He and C are 
fixed to their solar values. The remaining elements are combined into four groups for fitting: Mg; Si; S and Ar; and Ca, Fe, and Ni. 
Large systematic uncertainties in the ASCA SIS at low energies prevent us from determining the O abundance. 

A specific problem in the analysis of this sample, compared to the cluster samples of FAD and FDP, is the low signal-to-noise 
ratio for most of these groups. Given the spectral shape for 1 keV plasma, in most of the systems, no continuum is detected above 2 
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keV in the outer annuli. This limits the S abundance determination to only a few systems (NGC2300, NGC507 and NGC3268) and 
increases the uncertainties in the measurement of the Si abundance. Low signal-to-noise also makes the Mg abundance determination 
difficult, although temperature and Fe abundance determinations are robust. 

Heavy element abundance profiles are shown in Fig s .|7|— [To|. Like the temperature, elemental abundance profiles evolve with radius 
from stellar mass loss values to intragroup medium values, as indicated by an increase with radius in the gas-to-light ratio. Later types 
of the chemical enrichment are given preference in the (central) zones with low gas-to-light ratios, characterized by the dominance 
of SNe la in the Fe production. 

With the possible exception of the Mg abundance in NGC3268, the observed heavy element abundance profiles can be understood 
within the framework developed in FDP, FAD and Finoguenov et al. (2002). Groups are expected to have a lower abundance of 
SN II products as well as prevalence of SN la in the Fe production, related to low gas fraction in groups. The high gas temperature 
in NGC3268 causes the Mg measurements to be compromised (Mushotzky et al. 1996). This is further complicated by the strong 
temperature gradient in NGC3268, which, for the abundance determinations of light elements could not be sufficiently resolved in our 
data. Thus a confirmation of the Mg abundance profile requires observations with higher spatial resolution, such as XMM-Newton. 
Given the low X-ray luminosities of the systems in this study, determination of the central elemental abundance needs much higher 
spatial resolution in order to illuminate contributions from galactic point sources, a central AGN and a possible diffuse hard source 
(Finoguenov, Jones 2001). We therefore caution the reader about the limitations of our abundance results in the central spatial bins. 

Determination of the element abundance, characteristic of the elliptical galaxies has raised a problem when compared to optical 
data and raised controversy among the X-ray studies (Buote 2000, FDP). Element abundance gradients, reported here, are at least 
partly due to transition between the galaxy ISM and intragroup gas. While a number of pitfalls exist in determining of the element 
abundance for ellipticals from X-rays (Finoguenov et al. 2002), abundance determination, corresponding to the intragroup gas is 
robustly determined with ASCA data. 
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FIG. 7. — Derived iron abundances (in units of 4.68 1CT 5 for iron number abundances relative to H). The solid lines correspond to the best-fit Fe abundances derived 
from the ASCA data. The filled circles indicate the spatial binning used in the analysis. Dark and light shaded zones around the best fit curves denote the 68 and 90 
per cent confidence areas. 
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FIG. 8. — Derived silicon abundances. Data representation is similar to Fig. 
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FIG. 9. — Derived magnisium abundances. Data representation is similar to FigQ. 





FIG. 10. — Derived sulfur abundances. Data representation is similar to Fig.^J 
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Table 2 

ASCA SIS MEASUREMENTS OF TEMPERATURE, ENTROPY AND ELEMENT ABUNDANCE* 
Annulus(') kT, keV S.keVcm 2 Si/ Si® Fe/Fe G 





NGC5129 










7—15 


0.620 (0.58-0 


.65) 


280 ±19 


0.507(0.31-0.69) 


0.193 (0.17-0.21) 


3.2—7 


1.120(1.00-1 


.23) 


423 ±69 


0.465 (0.37-0.54) 


0.240 (0.15-0.33) 


1.5—3.2 


0.993 (0.94-1 


.03) 


162 ±17 


0.436 (0.32-0.54) 


0.299 (0.26-0.34) 


0—1.5 


0.914(0.87-0 


.94) 


77±5 


0.443 (0.25-0.63) 


0.423 (0.37-0.45) 




IC4296 










7—15 


0.931 (0.83-1 


.01) 


241 ±33 


0.318 (0.21-0.46) 


0.056 (0.03-0.09) 


3.2—7 


1.151 (1.03-1 


.34) 


308 ±83 


0.360 (0.20-0.58) 


0.323 (0.25-0.44) 


1.5—3.2 


1.110(1.02-1 


.32) 


177 ±5 


0.459(0.17-0.83) 


0.556 (0.42-0.74) 


0—1.5 


0.669 (0.61-0 


.72) 


56 ±14 


0.543 (0.15-1.06) 


0.687 (0.46-1.03) 




NGC4325 










4—14 


0.960 (0.83-1 


.16) 


435 ± 158 


0.718 (0.47-0.99) 


0.112(0.04-0.18) 


1.5 — 4 


1.244(1.12-1 


.40) 


166 ±29 


0.480 (0.33-0.65) 


0.244(0.16-0.37) 


0—1.5 


0.901 (0.88-0 


.93) 


39 ±3 


0.284 (0.18-0.43) 


0.296 (0.25-0.38) 




NGC7619 










5—10 


1.023 (0.98-1 


.07) 


200 ±16 


0.480 (0.36-0.63) 


0.209(0.18-0.25) 


2—5 


1.007 (0.71-1 


.28) 


336 ±327 


0.409 (0.26-0.55) 


0.212 (0.08-0.34) 


0—2 


0.806 (0.77-0 


.84) 


50±6 


0.339 (0.20-0.50) 


0.241 (0.20-0.28) 




NGC2300 










6—20 


1.215 (1.11-1 


.35) 


364 ±47 


0.390 (0.26-0.54) 


0.061 (0.03-0.10) 


2—6 


1.143 (0.92-1 


.38) 


215±73 


0.246 (0.10-0.40) 


0.101 (0.02-0.20) 


0—2 


0.819(0.75-0 


.88) 


45 ±10 


0.152(0.00-0.37) 


0.139 (0.09-0.21) 




NGC507 










4.3—8 


1.480(1.36-1 


.63) 


194 ±18 


0.373 (0.31-0.44) 


0.271 (0.24-0.30) 


2.3^.3 


1.293 (1.11-1 


.55) 


130 ±23 


0.563 (0.45-0.66) 


0.251 (0.21-0.28) 


1.2—2.3 


1.173 (1.07-1 


.33) 


62±7 


0.632 (0.57-0.70) 


0.258 (0.23-0.28) 


0—1.2 


0.772 (0.68-0 


.85) 


34±5 


0.599 (0.41-0.77) 


0.322 (0.28-0.36) 




NGC2563 










4.2—9 


1.441 (1.34-1 


.60) 


350 ±48 


0.670 (0.46-0.92) 


0.285 (0.20-0.39) 


2^1.2 


1.367 (1.22-1 


.58) 


231±53 


0.751 (0.59-0.97) 


0.291 (0.19-0.39) 


0—2 


0.983 (0.88-1 


.06) 


100 ±27 


0.854 (0.66-1.20) 


0.265 (0.20-0.44) 




NGC6329 










5.3—14 


1.485 (1.30-1 


.72) 


638 ±128 


0.496 (0.34-0.70) 


0.275 (0.21-0.37) 


2—5.3 


1.600(1.36-1 


.86) 


35 ±66 


0.327 (0.21-0.45) 


0.156 (0.10-0.24) 


0—2 


1.041 (0.98-1 


.09) 


82±9 


0.113(0.00-0.23) 


0.140(0.10-0.18) 




NGC3268 










24.2—35 


1.721 (1.60-1 


.77) 


445 ±22 


0.459 (0.42-0.49) 


0.100(0.08-0.14) 


16.7—24.2 2.195 (2.09-2 


.28) 


486 ±28 


0.493 (0.46-0.52) 


0.344 (0.31-0.39) 


11.6—16.7 2.588 (2.53-2 


.70) 


465 ±21 


0.496 (0.47-0.51) 


0.619 (0.57-0.65) 


8—11.6 


2.700 (2.60-2 


,79) 


493 ±39 


0.515 (0.48-0.56) 


0.886 (0.84-0.92) 


0—8 


2.625 (2.53-2 


.71) 


351 ±14 


0.533 (0.48-0.58) 


0.960 (0.91-1.00) 



' Definition of the solar abundance units is 3.55 and 4.68 X 10 5 for the number abundance of Si and Fe relative to H, respectively. Errors are given at 68% 
confidence level (see text for further details). MEKAL plasma code is used for spectral fitting. 



